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A cycle expansion technique for discrete sums of several PF operators, 
O" 

similar to the one used in standard classical dynamical zeta-function formal- 
ism is constructed. It is shown that the corresponding expansion coefficients 

\f\ ' show an interesting universal behavior, which illustrates the details of the 

(N' 

interference between the particlar mappings entering the sum. 

q' 
u 

pH ■ I. GENERAL INTRODUCTION 

Using the formalism of dynamical zeta functions one can compute global averages of the 
observables, associated with a fully chaotic dynamical system 

O ■ ii = Fi(x). (1) 

O' 

In order to find the average for an observable $(x), one defines a Perron-Frobenius operator, 

^ L^y) = 5{y,-f,{x))e^ x \ (2) 

> 

where f^(x) is a flow associated with the system (ffl), and f3 is a parameter. The time average 



ctf 



of $(x) can be found as 

($> = — lnz min (/3) l^o, (3) 

where z m in(/5) is the smallest root of the eigenvalue equation 

det (1 - zL) = 0. (4) 

For fully chaotic (axiom A) systems the determinant (Q) is an entire analytic function of 



z. 



which can be written [l], 2, 3] in an infinite product form as 



det (i - zl) =nn i- IAIAm - nn (i - ^^%) ■ ^ 

m p \ \ P\ p / m p 

Here t p = |A | Am is the weight of the prime periodic orbit of the system (P indexed by p, 
A p is it's Lyapunov exponent, A p is a short notation for the amplitude, and m enumerates 
repetitions of periodic orbits. The product over the prime periodic orbits only, 

Q\z) = II C 1 - z nv ^t p ) , (6) 

p 

is known || as the dynamical zeta function associated with the system ([!]), which is a 
meromorphic function in a certain domain of the complex variable z |], 2, 3]. 

An important feature of the theory based on the dynamical zeta function considerations, 
is the possibility to determine the coefficients of the power expansion of this function, 

C" 1 = 1 - tfz - c 2 z 2 - c 3 z 3 - ..., (7) 

using very effective and insightful "cycle expansion" technique ||]. As an example, in the 
case of a system whose orbits can described by a simple binary code, one can write the 
following cycle expansion of the zeta function C^ 1 - 



■vm 



C 1 = (1 - Zto) (1 - Zh) (l - Z 2 t 01 ) (l - 2 3 t 00 l) (l " Z 3 t 011/ 

= 1 — (ti+to) z — (toi — toti) z — ((tooi — ^o^oi) + (ton — titoi)) z - 
It is easy to observe that in the "curvature coefficients" , which are in this case 

C2 = toi — toti, 

c 3 = (^ooi — ^o^oi) + (^011 — ^l^Ol) ? 



(8) 



the contribution from long orbits is mimicked by the combinations of the short ones, which 
contribute to the sum with the opposite sign. As a result, the magnitude of the coefficients q 
of this expansion rapidly decrease. This allows to estimate very effectively the asymptotics 
of the coefficients of the expansion j|, 1], and to prove the analyticity properties stated 
above. 



II. GENERALIZATIONS 

The basis of the PF operators theory was developed by Grothendieck in [H] . His ideas 
were generalized recently by Ruelle || and Kitaev 0, who proposed to consider formal sums 
or integrals of PF operators L w , depending on a certain (discrete or continuous) parameter 

uj: 

Mx,y)= (i(dw)L u (x-f u (y)). (9) 



Here fj,(du) is some appropriate measure ||. One of the most important requirements 
imposed on such a sum or an integral, is that the corresponding dynamical systems have to 
be fully chaotic for all the values of u € Q. 

As shown in ||, under certain natural requirements imposed on the dynamical mapping 
function f(x) and the amplitude, the corresponding £ functions and the Fredholm determi- 
nants are analytic functions in a certain domain of the complex plane. However, the size 
of this domain is smaller than the one in the case of a single PF operator. Even if the 
individual systems L^ (f w ) have an entire Fredholm determinant Z and the zeta function, 
the size of the analyticity domain of the averaged operator (L) can be finite. 

Physically, operators like (|8|) can be used to study chaotic systems influenced by noise. 
Also, as shown recently in ||, the quantum mechanical Green's function of certain vector 
quantum particles can be presented as a sum similar to ©. Recent series of publications 
offered some effective methods for dealing with a system perturbed by Gaussian noise, 
including that of smooth conjugations |J and local matrix representation |10 |. However 



within these techniques the connection with the cycle expansion is not transparent. In this 
paper it is attempted to construct a cycle expansion technique for discrete versions of such 
sums, for which an analogue with cycle expansion technique is easy to establish. 



III. DOUBLE PF OPERATOR 

The simplest "generalized PF operator" is a formal sum of two (noncommuting) PF 
operators (J|), 

L®(x,y) = AS (y, - fl(x)) + A 2 5 (y, - f 2fM (x)) , (10) 

Here A\ and A 2 are certain multiplicative amplitudes. 

Following the standard procedure |l], 2], in order to describe the asymptotics of the 
evolution determined by the operator L^ 2 \x,y), we need to evaluate it's largest eigenvalue, 
or the smallest root of the equation 



oo z n 



det (l - zL^) = exp - £ -Tr (l&) = exp £ £ Tr (L h ....L in ) = 0, (11) 



n=\ U n=l \i\= n 



where ik = 1,2. If a sequence Li x ....L iri (i& = 1,2) is periodic, it can be written as a power 
of a shorter aperiodic string, 



L h ....L in = (L i /....Lj/J , (12) 

where n = s ■ q. Also, since only the trace of these operators is considered, the sequences 
obtained from one another by cyclic permutations give identical contributions, and therefore 
any length-s string contributes exactly s times. Hence 

J2Tr(L h ....L in ) = J2J2 s - Tr { L i'i-- Li 's) ' *i > -,in = l,2, (13) 

|i|=n g=l \i\=s 

and for the spectral determinant (pTDj) one has: 

oo oo yS-q 

log det (l - zL&) = E E E — ^ (^....L,,) 9 , (14) 

q=ls=l\i\=s $ 

where the sum is written in terms of the powers q of the aperiodic strings of the length s of 
the symbols L x and L 2 . 
Using the formula 



one gets the following expression for the trace of a generic term in the exponent: 



Tr{L ix ...L is ) q = j L h (x,yx) ...L in (y n -i,x) dxdy x ...dy n - X = J 6 (x - f h (y)) x ... (16) 



■•• x <5 (y ~ flS x )) dxdy 1 ...dy n - 1 = 6 (x - f h (y) o ... o f in {x)) = E 1 - 

J u 1-A 



(•) 

w,q 



(17) 



where w's are the fixed points of the map f^ (yi) o ... o fi n (x). The numbers s% and S2 
show how many times the operators L% and L 2 appear in the string L il ....L is . Obviously 
si + s 2 = s. Here for simplicity we put A\ — A 2 = 1, to avoid unnecessary complications 
in formulas. The dependence on the amplitudes is easy to reconstruct at the final step, so 
they will be ignored until the very last section. 

Hence, there is an infinite number of nonequivalent operators L^ ....L^ corresponding to 
each aperiodic string, and we are looking at their fixed points. Proceeding as usual, we can 
expand the denominator: 

1 ^ ^ 1 



1-A 



w,q 



EE 
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(s)m ' 
w,q 



(18) 



A V . 
w m=0 \-t'-w,q 

Here (s) denotes any of the aperiodic sequences of the length s, and A^ represents the 
expansion factor of a prime orbit indexed by "p" of a composite operator Lj 1 ....L, s . Substi- 
tuting this expansion series into the exponent yields the following Fredholm determinant: 

oo oo 



det(l-zL)=exp-££££- 

(s) V m=0r=l T 

Correspondingly, the product 



A 



W 



A 



(s)m 
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m=0 ( s ) p 
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{s) 



A ( 8 )m 



(19) 
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(20) 



can be considered as a generalization of the zeta function of the single-PF operator (^). 

This object is quite different from it's one-mapping counterpart. While the single- 
mapping zeta function tffty is defined on the periodic orbits of the system (|1|), the product 
(~ l is taken over the fixed points of the (aperiodic) combinations f# o ... o fy , which can be 
thought of as "interference" terms. Expanding the (s) product in (|20|) , one has 



c 1 = n i 



2\ n P 



A« 



A P ; 



(l)m 



1 



(z 2 ) 

A (2) A? )m 
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Ai 3) Af m 



(21) 



In order to avoid unnecessary complication of formulas, the repetition index m will be 
suppressed in the sequel. The parentheses (1), (2), (3) and so on imply the product over 
all the possible aperiodic combinations L^-.-.L^, i = 1,2, of lengths 1, 2, 3 correspondingly. 
For instance, the first bracket consists of the simple product of two (^-functions (two possible 
length- 1 strings), 



n (i - z n nu) = n (i - z n %) n f 1 - **>* 



(22) 



where t l = , j — . This is just the product of the individual (^-functions, 

\ p | p 



C 1 (z | U) = n (i - z n nl) = l - 1) 



J r 2 1 _3 

-2 Z C 3 X 



f-/^ c 2 z • c 3 z ..., 



c 1 (* I l 2 ) = n (i - *"*$ = i - 1 



! _ 2 2 _ 2 3 

c 2- Co Z Co -2- ... . 



The term t\ and t} represent the fundamental cycles, and the coefficients c\ and c 2 denote 
the curvature corrections for the first and the second mapping correspondingly. 

The only length-two string L\L 2 (the first "interference" correction) produces it's own 
zeta function, which contributes the term 



C 1 (z I L X L 2 ) = J] (l - z 2n Hf] = 1 - t 12 z 2 - cl 2 z 4 



4V... 



(23) 



to the product. Here t l ? , c 2 2 , C3 2 , etc., are the fundamental cycles and the curvature correc- 
tions for the mapping f\ o f 2 (x). 

The two length-three strings L 1 L 1 L 2 and L\L 2 L 2 contribute 



C 1 (z I l^l,) = n (1 - z 3n nf 2 ) = 1 



zHf 2 



c. 



112^6 ^112^9 



(24) 



and 



c 1 (« 1 ^^L^ = n (1 



n p .122 
6 p 



1 - zHf 2 



122 fi „122 9 



Co Z 



C 3 Z 



(25) 



correspondingly, and so on. The top index shows the aperiodic strings of operators L\ and 
L 2 , while the lower index enumerates the prime orbits of each operator, tf represents 

6 



X 



2\ 2 



the fundamental cycle of a given (s) operator, L il ....L is) and the coefficients c£> give the 
corresponding curvature corrections. 

It is important to notice, that the symbolic dynamics of the individual mappings, used 
to obtain the expressions for C" 1 (z \ Li), C _1 (z \ LiL 2 ), C _1 (z \ L\L\L<i), C,~ x (z \ L1L2L2), 
etc., do not need to be specified. In fact, one might use different symbolic dynamics in order 
to obtain the power expansion for these functions. After the necessary expansion series are 
at hand, the details of the symbolic dynamics of the corresponding mappings are irrelevant. 
The further analysis is made entirely in terms of the tf s and c£>s. 

Let us consider the power expansion of all the zeta-products up to the third power of z. 
The expansion for the two-PF operator will be then: 

(- l {z) = (l - t)z - c\z 2 - c l 3 z 3 - O (V)) (l - t)z - cjz 2 - cjz 3 - O (z 4 )\ 
(l _ z Hf _ o (z 4 )) (l - z 3 tf 2 - O (V)) (l - zHf 2 - O (z G 

= i-(t} + tj)z-((tf-t}tj) + 4 + 4 

[tf 2 -t 2 ftf)-(tr-ty f 2 )-( t } + t 2 f )(ci + ci) 

+t)c\ + t)c\ + c\ + 4] z 3 + (/) , (26) 

and so on, where tVs and d n are the fundamental contributions and the curvature corrections 
of the individual mappings. 

Comparing the expressions (^B|) and (§]||), one can see that the structure of the coef- 
ficients tf and c n is repeated on a "new level". The new "fundamental" contribution is 
now 

Tf = t} + tJ, (27) 

and the following terms have a structure similar to the corresponding "curvature corrections" 
in the case of a system with complete binary dynamics @ . For example the term 

C 2 = tf - t}t}, (28) 

is formally identical to the c-i of (ffi), the term C3 also contains terms 



{tr-t)tf)+{t)tf-tr) (29) 

of manifestly "curvature correction" type. 

However, starting from C3 there start to appear some additional "curvature combina- 
tions", such as 

(t) + t))(c\ + 4)-t)cl-t)ci (30) 

Also, every term contains the sum of all the standard curvature coefficients c n of the corre- 
sponding degree, associated with the "pure" maps fx(x) and fiix). 

In short notations, which are going to be used below, the expansion fl26|) can be written 
as: 

C2\ / 2 2 2 2\ 

8=1 / \ 8=1 8=1 8=1 8=1 / 

(31) 

It should be emphasized, that the coefficients Tf and C n are principally different from their 
usual, single PF counterparts tf and c n . They are constructed using the contributions of 
different mappings, (not only the fi(x) and ^(a;) but also from all their aperiodic superpo- 
sitions), and therefore produce some "interference" effect. 

IV. TRIPLE PF OPERATOR 

One can easily write down the explicit form of the power expansion for the triple PF 
operator 

L^(x,y) = A x 8 (y, - fax)) + A 2 5 (y„ - fax)) + A 3 5 (y„ - fax)) . (32) 

As in the double PF operator case, the expansion is made up to the third order: 

C3- 1 = n (1 - t)z - 4z 2 - c 3 , 3 - o(z 4 )) n (1 - z hj - o^)) n (1 - Z Hf - o^ 

i=l ij ijk 



l-z(t}+t} + t}) - z* ((tf - tf}) + (tf - tf}) + (tf - tf}) +cl + cj + c 



+ z 3 



122 _ ,2 +12 
/ l f l f 



i /.113 +1+13\ | ?+!33 +3+13\ , .213 , +132 +3+12 .2.13 +1+23 , +123 , +1+2. 

+ Vf ~ t f t f ) + V / ~ / / ) + f + f ~ f f ~ f f ~ J f + / + / / 



-t) (4 + 4) - 1) (c\ + cf) - 1) (c\ + 4) + (c\ + 4 + 



+ 0(A (33) 



where tf, c n are the fundamental contributions and the curvature corrections of the corre- 
sponding mappings. 

It is easy to observe that the above expansion is very similar to the one obtained for a 
system with a complete ternary symbolic dynamics, for which the expansion coefficients are: 

4 3) = *s + h + t x 

cf = («13 - M 3 ) + (*23 " hh) + (tl2 " ht 



1«-2J 



,(3) 



C3 — (^233 — M23) + (^223 — M23) + (^133 — M13) + (^113 — M13) + (^122 — M12) 
+ (^112 — M12) + ^213 + ^123 — Ml2 — Ml3 — M23 + ^1*2^3 

(3) (3) (3) 

Using the notations T) , C2 , C3 for the corresponding combinations of prime contri- 
butions in (PEP, we get for the triple operator zeta functions: 

C3- 1 = 1 - zT f - z* (cf + E 4) - ^ 3 (c? - E t) E 4 + E 44 + E 4) + o(A 

\ i=l/ \ i=l i=l i=l »=1 / 

(34) 
Comparing this to the expansion of the double PF operator, 

C2- 1 = 1 - *Tf - , 2 (cf + E 4) - - 3 f cf ] - E 4 E 4 + E 44 + E 4) + o(A 

V i=l/ V t=l i=l i=l i=l / 

(35) 

it is easy to observe that the curvature coefficients have the same structure, depending 
trivially on the number of the PF operators in the sum elements. The main difference is 
contained in the terms C 2 , C3 as opposed to C^ , C3 , because they are constructed as 
the curvature terms of systems with binary and ternary symbolic dynamics correspondingly. 



V. tV-SUMS AND THE CONTINUOUS LIMIT 

Previous analysis can be conducted in the exact same way for any iV -operator sum 
(N>2), 

TV 

L^(x,y) = Y,S(y,-fi,(x)). (36) 

*=i 

For the Fredholm determinant det (1 — zL) one has: 

det(i-,L)=ex P -E-^(^ n ) =nnn i- na J - ( 37 ) 



77 

ra=l / m s p 



A 



A 



•p ilp 



Based on the previous examples and using the induction method, one can prove that the 
power expansion of the zeta function of an iV-term PF operator is 



G 1 (,) = i-,Tf)-, 2 (c7f) + E4) 



/ TV TV TV TV \ 

-z s (ci N) + ^f f J2ci -E44 + E4 + 0(z% (38) 

\ 8=1 8=1 8=1 8=1 / 

where 

TV 

if = £4 (39) 



and C\ , C3 , etc., have the structure identical to the curvature correction coefficients of 
a system with a complete iV-ary symbolic dynamics. The extra terms which appear in the 
expansion (p8|), 

TV 



(40) 



(41) 

can be called "direct curvature contributions" as opposed to the "interference curvatures" 

C( N ). 
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I — 1 

TV TV 
i=\ i=\ 


TV TV 

-E44 + E4 

8=1 i=\ 



The direct curvature terms fl40|) depend trivially on the number of PF operators in- 
cluded into the sum fl3"T|). Among the interference curvature terms, only the fundamental 
contribution term Tf is just the direct sum of the individual fundamental contributions. 

It is important to emphasize, that the complete curvature coefficients have a universal 
form, 

Ci N) = C^ + c^\ (42) 

which provides an algorithm for evaluating the zeta-function expansion coefficients of the 
multiple PF operators to a given degree n. Moreover, the direct curvature coefficients depend 
trivially on the order of the sum (^). This allows to make certain statements about the 
behavior of the continuous sums, 

L(x, y)= J n (duo) A J (y„ - /* (x, uj) , (43) 

from the point of view of the standard cycle expansion technique. Here the amplitudes A^ 
and the mapping functions /* (x,u) depend continuously on the parameter u. 

After some elementary considerations it is easy to reconstruct the explicit amplitude 
dependence in the expansions of the MPF operators' zeta functions, 

/ 7 s-n p A(s)n p \ 



A P 



A P 



(s) v 

where A^ s ' = A S i ...A S 2 N . The product over the (s) includes all the possible strings in which 
the symbols L lr .., L N appear Si,..., s 2 times correspondingly, for all s\ + ... + s^ = s. The 
expansion of (|44]) yields: 



N ( N \ 

8=1 V 1=1 / 

/ N N N N \ 

-z* a®c 3 + Y, ah - E A lt ) y: 4 2 4 + E A 3 44 + o(z% 

\ i=l i=l i=l »=1 / 

Without loss of generality, one can assume that there exists such a set of discrete pa- 
rameters u>x, ...,un, that the curvature corrections of the individual maps can be written 
as: 

11 



t) = tf (Ui) , 



4 = c n (ui) . (45) 

In case if the amplitudes A^ and the mapping functions /* (x, to) in fl44| ) depend continuously 



on the parameter u, the coefficients (%3 ) also depend continuously on w, at least for a certain 
range of parameters. Therefore, the fundamental contribution coefficient will produce in the 
limit N — > oo: 

AT . 

if ] = £ A4 - 2> = / t f ( U )ii (dco) . (46) 

4 = 1 

The direct curvature contributions in the expansion fl40| ) also can be written in the integral 
form: 

N N 



ci N) = Y: Ajcl = Y: A*c 2 (w,) - £ 2 = / c 2 (w) /x (d 2 u;) , 

«=1 «=1 

c 3 = / tf(w)(i, (dco) I c 2 (io)n (d 2 ujj - / ^02/2 (d 3 u;J + / c 3 (u)fi (d 3 oo) , 



(47) 

On the other hand, the interference curvature coefficients C^> the limit iV — > oo are 
considerably more complicated, since the elements contained in the higher order interference 
curvatures C^ N \ n > 3, are labeled by the aperiodic strings of complexity N . 

However, the N — > oo limit for the interference coefficients exist P,|TO|, and therefore 
zeta function C _1 has the following expansion form: 

Q 1 ^ = l-zf t f (u)fi {duo) - z 2 (C» + Jc 2 (oo) fl (d 2 oofj 
-z 3 I Cf + J t f {uo)n (duo) J c 2 (lo)io (d 2 uo) + J t f c 2 (uo)ii (d 3 uo) 

-Jc 3 (oo)fi(d 3 oo))+O(z 4 ). (48) 



This expansion is analogous to the cumulant expansion of det (1 — zL) obtained in jjlO . 
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VI. CONCLUSION 



From practical as well as from theoretical point of view, it is certainly more convenient to 



treat the continuous sum case using the methods developed in P,|10|, because these methods 
effectively go around some convoluted calculations, which are necessary to evaluate the 
expansion coefficients (pH). However, for discrete versions of such "noisy" PF operators, the 
multiple sums (^), these methods can not be applied directly. Meanwhile, as it was shown 
above, such systems allow a direct treatment which is analogous to the standard cycle 



expansion technique. The detailed form of the corresponding expansion coefficients (ID) 



reveals a remarkable structural simplicity and shows certain a familiar features characteristic 
for the single-operator expansion coefficients. 

In practice, even in case of the a single PF operator, one can usually obtain only the first 
several terms of the cycle expansion series. In case of multiple PF operators, the explicit 
form of the generalized curvatures ( [40|) conveniently illustrates the details of the interference 
between the different maps constituting a multiple PF operator, using the language of the 
periodic cycles arising from the maps themselves as well as their compositions. This provides 
a nontrivial extension of the periodic orbit theory to a more complicated dynamical evolution 
which involves a few "evolution channels". 
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